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We apply the density-matrix renormalization group technique to magnetic molecules in order to 
evaluate therjow-lying energy spectrum. In particular, we investigate the giant Keplerate molecule 
{Mo7 2 Fe 3 o},El where 30 Fe 3+ ions (spins 5/2) occupy the sites of an icosidodecahedron and interact 
via nearest-neighbor antiferromagnetic Heisenberg exchange. 

The aim of our investigation is to verify the applicability and feasibility of DMRG calculations 
for complex magnetic molecules. To this end we first use a fictitious molecule with the same 
structure as {Mo72Fe3o} but with spins 1/2 as a test system. Here wpiinvestigate the accuracy of our 
DMRG implementation in comparison to numerically exact resultsu Then we apply the algorithm 
to {Mo72Fe3o} and calculate an approximation of the lowest energy levels in the subspaces of total 
magnetic quantum number. The results prove the existence of a lowest rotational band, which was 
predicted in Ref. |j| 
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I. INTRODUCTION 

The rapid progress in polyoxometalate chemistry gen- 
erates larger and larger magnetic molecules. The most 
prominent example of recent times is the molecular 
magnet {Mo72Fe3o}, where 30 Fe 3+ paramagnetic ions 
(spin 5/2) occupy the sites of an icosidodecahedron, see 
Fig. |l|, and interact via isotropic nearest-neighbor anti- 
ferromagnetic Heisenberg exchange Whereas the sta- 
tistical and dynamical properties of smaller molecules 
can be evaluated by numerically exact diagonalization, 
the huge dimension of the Hilbert space of magnetic 
macromolecules prohibits such attempts. In the case of 
{Mo72Fe3o} this dimension amounts to 6 30 , which is be- 
yond the power of any computer. 




FIG. 1: Three-dimensional model of {Mo72Fe3o}: the vertices 
mark the sites of the Fe 3+ ions (spin 5/2), and the lines denote 
nearest-neighbor interactions. 



Many attempts have been undertaken in order to ap- 
proximate the energy eigenvalue spectrum of large mag- 
netic systems. Among these the Density Matrix Renor- 
malization Group technique (DMRG)H is one of the most 
powerful, at least in the field of one-dimensional spin sys- 
tems, where for instance questions concerning Haldane's 



conjectureBLJ could be answered with great accuracyJ3&E2l 
It is as well a powerfuL-tool for studying T = quan- 
tum phase transitions]^ as it delivers accurate results 
for ground states. The DMRG method was also ap- 
plied to spin rings - "ferric wheels'?— ■ which are quasi 
one-dimensional magnetic moleculesJ13 Here the aim was 
to evaluate low-lying magnetic levels and the related 
low-temperature spin dynamics in order to understand 
macroscopic quantum coherent phenomena. 

Another method of approximating the energy spec- 
trum is stimulated by the observation that in many 
Heisenberg spin systems the low-lying energy levels 
E m in(S) form a rotational bandJ13 i.e. they depend 
approximately quadratically on the total spin quantum 
number S. Experimentally this property has , 
scribed as "following the Lande interval rule" olljlif 
For spin ring systems the low-energy spectrum consists 
of a sequence of rotational bands, which allows to address 
questions of spin tunneling and other, transitions without 
diagonalizing the full Hamiltonian.EEl 

The purpose of this article is twofold. We show 
that DMRG can be used in order to approximate the 
low-lying energy levels of magnetic macromolecules like 
{Mo72Fe3o}, and we prove numerically that the lowest 
levels as a function of total spin S form indeed a ro- 
tational band. The latter observation strengthens the 
predictions made in Ref. 0. 



II. DMRG TECHNIQUE 

The DMRG technique!! became one of the standard 
numerical methods for quantum lattice calculations in 
recent yearsJij Its basic idea is the reduction of Hilbert 
space while focusing on the accuracy of a target state. 
For this purpose the system is divided into subunits - 
blocks - which are represented by reduced sets of ba- 
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sis states. The size m of the truncated block Hilbert 
space is a major input parameter of the method and to 
a large extent determines its accuracy. The block ba- 
sis states are derived from a twice as large system - su- 
perblock - by first diagonalizing the Hamiltonian on the 
superblock, then building a reduced density matrix from 
the superblock ground state, and finally diagonalizing the 
reduced density matrix. 

DMRG is best suited for chain-like structures. Many 
accurate results have been achieved by applying BjMRG 
to various (quasi-)one-dimensional systems .l3tJ't£l The 
best results were found for the limit of infinite chains with 
open boundary conditions. It is commonly accepted that 
DMRG reaches maximum accuracy when it is applied to 
systems with a small number of interactions between the 
blocks, e.g. systems with only nearest-neighbor interac- 
tion and open boundary conditionals. 

In order to apply DMRG calculations to two- 
dimensional systems a mapping onto a one-dimensional 
structure was proposed.Ea We adopt this idea and de- 
rive a simple DMRG algorithm hp-two-dimensional spin 
systems in the Heisenberg model.EHI Since the spin array 
consists of a countable number of spins, any arbitrary 
numbering is already a mapping onto a one-dimensional 
structure. However, even if the original system had only 
nearest-neighbor exchange, the new one-dimensional sys- 
tem has many long-range interactions depending on the 
way the spins are enumerated, see e.g. Fig.|. Therefore, 
a numbering which minimizes long range interactions is 
preferable. 




FIG. 2: One-dimensional projection of the icosidodecahedron: 
the lines represent interactions. 

The Hamiltonian of the Heisenberg model, which is 
appropriate for the investigated magnetic molecules, can 
be written as 

% = - zJ J » & ■ & = ~ 2 J y ■ &■ ' w 

ijij i>j 

where Jij is the interaction matrix and S, are the spin 
operators at the sites i. 



1 2 • ■ 


p 


p+1 


• • N-1 N 


OO" 


• • o o • 


• OO" 


■ • o o 



FIG. 3: Block setup for second part of DMRG algorithm 
("sweeping"): The whole system of N spins constitutes the 
superblock. The spins {1, 2, . . . ,p} belong to the left block, 
the other spins ... , N} to the right block. 



We use a block algorithm similar to the setup in 
White's original paper.El Two blocks are connected via 
two single spin sites, these four parts form the superblock, 
see Fig. 0. The Hamiltonian for the superblock can be 
written as 

p N 

H SB = H l + Jf + (-2) £ £ JtjSi.Sj , (2) 
i=l j= P +l 

where H and H r represent the Hamiltonians for the 
left and right block including the respective single spin. 
These Hamiltonians include the interactions inside the re- 
spective blocks, therefore the third term in (|J) describes 
the interactions of the spins belonging to the left block 
with those of the right block. 

In our implementation we have to keep track of the op- 
erators Sf and Sf (in matrix representation) for all sites 
i. In the case of a system with pure nearest- neighbor in- 
teraction one would have to keep only the operators for 
sites at the borders of the blocks. Because of the long- 
range interactions, the 2D-DMRG approach consumes 
more memory and the calculation of the Hamiltonian 
takes more time. 

The algorithm consists of two steps. During the first 
step the superblock grows with each iteration by two sites 
until the final length of the system is reached. The second 
step is an implementation of White's sweep-algorithm.El 
While the superblock is kept at its maximum length, in 
each iteration the left block grows by one site whereas 
the right block is shortened by one site. When the right 
block reaches the size of two sites, the direction of the 
sweep is turned, and the right block grows in the next 
iterations. One performs a number of sweeps until the 
desired property, in our case the ground state energy, 
converges. 

The Hamiltonian is invariant under rotations in spin 
space. Therefore the total magnetic quantum number M 
is a good quantum number and we can perform our cal- 
culation in each orthogonal subspace H (M) separately. 

III. ACCURACY OF THE METHOD 

Since it is difficult to predict the accuracy of a DMRG 
calculation, we apply our implementation to an ex- 
actly diagonalizable system first. The most realistic 
test system for the use of DMRG for {Mo72Fe3o} is the 
icosidodecahedron with spins s — 1/2. This fictitious 
molecule, which possibly may be synthesized with vana- 
dium ions instead of iron ions, has the same structure 
as {Mo72Fe3o}, but the smaller spin quantum number 
reduces the dimension of the Hilbert space significantly. 
Therefore a numerically exact diagoiialization is possible 
and was carried out by J. RichtcrBEil We use these re- 
sults to analyze the principle feasibility and the accuracy 
of the method. 

The DMRG calculations were implemented using the 
enumeration of the spin sites as shown in Figs. and H. 
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high accuracy of one-dimensional calculations (often bet- 
ter than 10 -6 ), the result demonstrates that DMRG is 
applicable to finite 2D spin systems. 

Our results were obtained keeping m = 60 states per 
block. The number of states could easily be increased 
for the s — 1/2 case, but we wanted to have a prediction 
for {Mo72Fe3o}, where m is limited by the available com- 
puter resources because of the much larger spin s = 5/2. 
The number of sweeps ranged from 5 to 20 depending on 
how quickly the algorithm ran into oscillations with no 
further improvement of accuracy. This oscillatory behav- 
ior of the sweep algorithm is described in Ref. [L9. 



IV. ROTATIONAL BAND IN {Mo 7 2Fe 30 } 



FIG. 4: Two-dimensional projection of the icosidodecahedron, 
the site numbers are those used in our DMRG algorithm. 



We choose this enumeration because it minimizes the 
average interaction length between two sites. The inter- 
action length on the chain is measured as the distance 
\j — i\ of the two interacting sites i and j. The DMRG 
method favors systems with a minimal number of inter- 
actions between the blocks. Therefore, a short average 
interaction length helps to reduce the number of inter- 
block interactions. 

With our enumeration we get an average length of 3 
between two interacting sites. However, the choice is not 
unique because of the various symmetries of the system. 
Our choice, Fig. ^, is rotationally symmetric with a five- 
fold symmetry corresponding to the five-fold symmetry 
of the central pentagon. The sites 1 to 6 form the unit 
cell. 



Q exact eigenvalues 
DMRG calculation 
fit to rotational band 



FIG. 5: Eigenvalues and lowest rotational band of the s 
icosidodecahedron. 
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In Fig. H our DMRG results (crosses) are compared to 
the energy eigenvalues (circles) deterroiped numerically 
by J. Richter with a Lanczos method.00 We find very 
good agreement of both sequences, with a maximal rel- 
ative error of about 1%. Although we don't achieve the 



Since the DMRG technique has proven applicable for 
the s = i case of the icosidodecahedron, we use our algo- 
rithm to approximate energy eigenvalues of the magnetic 
Keplerate molecule {Mo72Fe 30 }. 

In Ref. U it was predicted that the low-lying energy 
eigenstates of {Mo72Fe3o} form "rotational bands", i.e. 
the sequence of ground states energies of the sub-spaces 
Ti. (S) is expected to have a quadratic dependence on the 
total spin quantum number S. A spectrum with ro- 
tational bands usually arises in antiferromagnets if the 
spin system can be divided into sub-lattices. The most 
prominent example are bipartite rings or chains which 
consist of two sub-lattices with opposite sub-lattice mag- 
netization. In the case of {Mo72Fe3o} the spin system 
is decomposable into three sub-lattices with_sub-lattice 
spin quantum numbers Sa, Sb, and Sc-Effij Then the 
low-lying spectrum can be described by an approximate 
Hamilton operator 



H --J- 

-2 approx — J jy 
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(3) 



where S is the total spin operator and the others are sub- 
lattice spin operators. At least in the case of bipartite 
systems this approximation has turned out to be a very 
good one£3ii3 

The minimal energy eigenvalues of H approx as a func- 
tion of S form a rotational band by construction 
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(4) 



We use the DMRG method to approximate the lowest en- 
ergy eigenvalues of the full Hamiltonian (0) and compare 
them to those predicted by the rotational band hypoth- 
esis (0). In our calculation we obtain energy levels for 
the ground states of TL (M) sub-spaces. These states are 
equivalent to the ground states of the sub-spaces H (S) 
with S = M. The proof for this property rests on the 
monotonous increase of the sequence E (M) with M for 
< M < Ns. 

Fig. ^| shows our results and a fit to the lowest ro- 
tational band. We find a good agreement between our 
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FIG. 6: DMRG eigenvalues and lowest rotational band of the 
s = 5/2 icosidodecahedron. 

DMRG data and the predicted quadratic dependence. 
From the fit of our data we obtain D = 6.17 and 7 = 1.05. 
These values are very close to the values D = 6.23 and 
7 = 1.07 given in Ref. ^, which were inferred from mag- 
netization measurements. 



molecule of 30 spins s = 1/2 and icosidodecahedral struc- 
ture. Therefore, we assume that our numerical approxi- 
mation of low- lying energy levels for {Mo72Fe3o} is reli- 
able. Thus, we have obtained good confidence that the 
prediction of a lowest rotational band made in Ref. ^ 
is justified. The lowest band of {Mo72Fe3o} indeed has 
a parabolic dependence on S. It remains the task of 
forthcoming investigations whether also the higher ly- 
ing rotational bands appear with the same distinctness 
or whether they are scattered due to the strong frustra- 
tion effects. In any case such calculations demand much 
higher precision and thus numerical effort. 

The present calculations were carried out keeping m — 
60 block states, which means that the calculation time for 
one ground state is about a day on a standard PC, and 
there are 76 Tt (M) ground states (M — 0, 1, ... , 75) in 
{Mo72Fe3o}. For a more accurate calculation one would 
have to use more powerful machines and more computer 
time both allowing higher values for m and therefore lead- 
ing to a better accuracy. 
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